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ABSTRACT 

The habitable zones of main sequence stars have traditionally been defined 
as the range of orbits that intercept the appropriate amount of stellar flux to 
permit surface water on a planet. Terrestrial exoplanets discovered to orbit M 
stars in these zones, which are close-in due to decreased stellar luminosity, may 
also undergo significant tidal heating. Tidal heating may span a wide range for 
terrestrial exoplanets and may significantly affect conditions near the surface. 
For example, if heating rates on an exoplanet are near or greater than that 
on Io (where tides drive volcanism that resurface the planet at least every 1 
Myr) and produce similar surface conditions, then the development of life seems 
unlikely. On the other hand, if the tidal heating rate is less than the minimum 
to initiate plate tectonics, then CO2 may not be recycled through subduction, 
leading to a runaway greenhouse that sterilizes the planet. These two cases 
represent potential boundaries to habitability and are presented along with the 
range of the traditional habitable zone for main sequence, low-mass stars. We 
propose a revised habitable zone that incorporates both stellar insolation and 
tidal heating. We apply these criteria to GJ 581 d and find that it is in the 
traditional habitable zone, but its tidal heating alone may be insufficient for 
plate tectonics. 

Subject headings: astrobiology — (stars:) planetary systems — stars: individual 
(GJ 581) — stars: low mass 



1. Introduction 

The discovery of extrasolar planets has made detecting and recognizing life-bearing 
planets outside the solar system a real possibility. Considerations of stellar radiation and 
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climate led to the definition of a "habitable zone" (HZ) as a region around a star in which a 
planet with an Earth-like atmosphere could support liquid water on its surface (Kasting et 
al. 1993; Selsis et al. 2007; Spiegel et al. 2009). Close to the star, the surface temperature is 
too hot for liquid water; far from the star it is too cold; but in between is a range of orbits 
for which liquid water is stable on the surface. This definition predicts that habitable orbits 
lie closer (~ 0.1 AU) to low-mass stars than to Sun-like stars because low-mass stars are less 
luminous. We call this range of orbits the "insolation habitable zone" (IHZ). 

Also of critical importance for determining habitability is the maximum mass that does 
not accrete a significant hydrogen-rich envelope. Numerous investigations have modeled this 
cut-off, and values range from 0.1 - 16 Earth masses (Me) (Pollack et al. 1996; Ikoma et 
al. 2001; Lissauer et al. 2009). For the expected conditions of most proto-planetary disks, 
these results suggest that 10 Me is a reasonable approximation for this critical mass, and 
hence, as is now standard, we assume that habitable planets must be less massive than this 
value. 

Transit and radial velocity surveys are most sensitive to planets in small orbits ( <0.5 
AU), and hence are unlikely to find terrestrial planets in the IHZs of solar-mass stars. How- 
ever, terrestrial-scale planets are readily detectable in the IHZs of low-mass stars (Mayor et 
al. 2009). Radial velocity surveys of M stars are complicated by the presence of stellar ac- 
tivity (West et al. 2004), which can induce surface fluctuations that erase detectable signals 
of terrestrial-scale planets. Nevertheless, the first plausibly terrestrial-sized exoplanet in the 
IHZ has just been discovered (GJ 581 d, > 7.1 M E ) (Mayor et al. 2009). Analogously, transit 
detections are hampered by starspots, which are more frequent and larger on M stars than 
on larger stars. Nonetheless, transit surveys may be the most effective method to discover 
terrestrial exoplanets due to their ability to break the mass-inclination degeneracy of radial 
velocity data, and measure the planetary radius directly. With combined transit and radial 
velocity data, the unambiguous detection of a terrestrial planet in the IHZ is feasible in the 
immediate future. Not surprisingly, surveys such as "MEarth" have begun to scan M stars 
for transits (Nutzman & Charbonneau 2008). 

However, planets within a few tenths of an AU from M stars are in a very different 
environment than anywhere in the IHZ of our Solar System: They may be hit with solar 
flares, bathed in X-rays due to increased stellar activity, and/or subjected to strong tidal 
effects. The first two are certainly concerns for the success of biology, although they were 
neglected in the original IHZ definition. Significant research has since explored these issues 
{e.g. Lammer et al. 2007; Segura et al. 2009; Zendejas et al. , in prep.). Here we focus on the 
third issue and show that tides may place important constraints on planetary habitability. 

For example, the recently-announced 2 Me planet GJ 581 e (Mayor et al. 2009) may 
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experience intense tidal heating. Mayor et al. showed that this planet's eccentricity occasion- 
ally reaches values of 0.1. Applying common models of tidal heating (e.g. Peale et al. 1979; 
Jackson et al. 2008a,b; Barnes et al. 2009), and assuming the planet to be terrestrial-like, GJ 
581 e could have 2 orders of magnitude more tidal heating than Jupiter's volcanic satellite 
Io! Although this planet is not in the IHZ, similar heating rates on planets in the IHZ are 
unlikely to develop life (Jackson et al. 2008a). 

However, only in some cases will the heating be greater than Io's; in others it may 
lead to geophysical processes (e.g. plate tectonics) that maintain long-term climate stability, 
increasing potential habitability (Williams et al. 1997). Jackson et al. (2008a) showed that 
tidal heating in the center of the IHZ could be substantial for a few illustrative cases. Here we 
identify the boundaries of this "tidal HZ" (THZ) for a range of stellar and planetary masses, 
M*, and M p , respectively, as well as semi-major axis a and orbital eccentricity e. This 
definition excludes radiogenic sources of heat, which dominate the heating on the Earth. 
Estimating radiogenic heating rates on exoplanets seems a daunting task as it depends 
sensitively on the composition of the planet (itself a result of many different formation 
pathways [Raymond et al. 2008]), the age, and the internal structure of the planet. Therefore, 
bear in mind that the heating rates we present here are in addition to any radiogenic heating. 

We define a new HZ to take into account tidal heat: Planetary habitability requires 
liquid water on the surface and enough internal heat to drive plate tectonics, but not so 
much as to cause Io-like volcanism. Combined with the other restriction, these additional 
requirements significantly reduce the range of detectable habitable environments around low- 
mass stars. In § 2 we describe the IHZ and THZ. In § 3 we combine these models to introduce 
the "insolation-tidal HZ" (ITHZ) and define it for the cases of 1 and 10 Me planets. In § 4 
we discuss the implications, and suggest directions for future research. 

2. Tidal and Atmospheric Models 

2.1. The Insolation Habitable Zone 

Kasting et al. (1993) suggested that planetary habitability required liquid water to be 
stable on a planetary surface. They made simple assumptions about the atmosphere, and 
the planetary orbit, (e.g. it was circular). Selsis et al. (2007) included the effects of cloud 
cover, and applied an updated version of the Kasting et al. model to the GJ 581 system 
(Udry et al. 2007). However, many of the known exoplanets are on eccentric orbits (Butler 
et al. 2006), implying rocky exoplanets will be also. Hence, the effect of eccentricity may be 
important when assessing habitability. Williams Sz Pollard (2002) pointed out that planets 
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on eccentric orbits tend to have surface temperatures that reflect the orbit-averaged stellar 
flux. 

Barnes et al. (2008) modified the Selsis et al. (2007) model to reflect the Williams & 
Pollard (2002) results in order to define an "eccentric HZ." Here we adopt the Barnes et 
al. definition of the IHZ, which requires a planet's semi-major axis to lie in the range U n < 
a < Lut, where l{ n and l ou t are the inner and outer edge of the IHZ, respectively. As in 
Barnes et al. , we assume all planets have 50% cloud cover (see Selsis et al. 2007). 



2.2. The Tidal Habitable Zone 

The exact mechanisms of tidal dissipation are poorly understood, but several quantita- 
tive models have been suggested (e.g. Goldreich & Soter 1966; Hut 1981; Levrard et al. 2006; 
Ferraz-Mello et al. 2008). A conventional model quantifies the tidal heating of a body as 

= 63 

4 Q' P y J 

where G is the gravitational constant, R p is the planetary radius, and Q' p is the "tidal 
dissipation function" which encapsulates the physical response of the body to tides, including 
the Love number (Peale et al. 1979; Jackson 2008b). Note that this model may break down 
for large e. With the exception of Q' p , all these quantities can be measured for exoplanets. 
However, in lieu of actual mass and radius measurements we use the scaling relationship 
R p oc M°- 27 (Sotin et al. 2007). 

In order to assess the surface effects of tidal heating on a potential biosphere, consider 
the heating flux, h = H/AitRl, through the planetary surface. On Io, h = 2 W m~ 2 (from 
tidal heating) (McEwen et al. 2004), resulting in intense global volcanism and a lithosphere 
recycling timescale of 142 - 3.6 x 10 5 years (Blaney et al. 1995; McEwen et al. 2004). Such 
rapid resurfacing on an exoplanet probably precludes the development of a biosphere, hence 
we will assume that heating rates larger than h max = 2 W m~ 2 result in uninhabitable 
planets. 

Internal heating can also drive plate tectonics. This phenomenon may enable planetary 
habitability because it drives the carbonate-silicate cycle, stabilizing atmospheric temper- 
atures and CO2 levels on timescales of hundreds of millions of years. Although the pro- 
cesses that drive plate tectonics are not fully understood on the Earth (Walker et al. 1981; 
Regenauer-Lieb et al. 2001) an adequate internal heat source is essential. Theoretical studies 
of Martian geophysics suggest tectonic activity ceased when the radiogenic heating flux h ra d 
dropped below 0.04 W m -2 (Williams et al. 1997). We will therefore assume the planets 
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have Earth-like weathering and that heating rates smaller than h m i n = 0.04 W m 2 result 
in uninhabitable planets. 

On Earth, the heating comes from the radioactive decay of U and K and produces a 
heat flux of /i m <2© = 0.08 W m~ 2 (Davies 1999). In the case of exoplanets, the radiogenic 
heating rates are unknown, but they might well be inadequate. Planets significantly smaller 
than the Earth, which may be common around M stars (Raymond et al. 2007a), may be 
especially susceptible to insufficient heating (the radiogenic heat flux scales as the ratio of 
volume to area). In many cases, tidal heating probably dominates (Jackson et al. 2008a,b), 
and may be directly estimated from the mass and orbit of the planet. Such a calculation 
helps assess a planet's internal heating and accompanying geophysical processes. 

The above examples depend critically on the internal structure and composition of the 
body in question. Without this information, objects in the solar system serve as the best 
guide for hypothesizing the internal structure and dynamics of exoplanets. Hence, while 
explicitly acknowledging the wide range of planets likely to be discovered, we consider the 
above definitions of h min and h max to be the limits of habitability. 

Therefore we define the THZ to be the region around a star for which h m i n < h < h max 
(where h only represents the tidal heating), in this case, enough heat for plate tectonics, but 
not so much as to result in extremely rapid resurfacing. As with the IHZ, this definition 
makes critical assumptions about the planetary properties, specifically the planetary Q' p 
value. We assume Q' p = 500, consistent with the values of most terrestrial-like bodies in the 
Solar System (Dickey et al. 1994; Yoder 1995; Mardling & Lin 2004). 



3. The Insolation-Tidal Habitable Zone 

Habitable conditions may be present on a planet if it lies in both the IHZ and the THZ. 
In other words, U n < a < l out and h min < h < h max . When both these requirements are met, 
a planet is in the "insolation-tidal HZ" (ITHZ). In this section we present examples of how 
this refinement modifies our view of habitability of planets orbiting main sequence stars. 

Figure [1] shows the ITHZ as a function of M* and a for a 10 Me planet, the most 
detectable terrestrial-like exoplanet. The blue shading represents the IHZ (see § 2.1); the 
yellow shading represents the THZ (see § 2.2). For the latter, two cases are shown: e = 
0.01 (left strip), and e = 0.5 (right strip). The intersections of the IHZ and the THZs are 
shaded green and are two ITHZs. Apparently by chance, the two habitability requirements 
significantly overlap in these and many other plausible cases, but the range of habitable 
orbits is smaller when considering both tidal and insolation effects. 
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In order to assess the likelihood of finding bodies in the ITHZ, we show transit detection 
probabilities as the red lines in Fig. [TJ The horizontal ones are geometric transit probabilities 
(calculated from Barnes [2007]), and the diagonal ones represent the percentage of starlight 
blocked during transit (the "transit depth", the ratio of the cross-sectional areas of the 
planet and star). Geometric probabilities increase with decreasing a, but transit depths are 
independent of a. Transit depths of ~ 0.5% have been detected from the ground (e.g. Winn 
et al. 2008). Transiting terrestrial-scale planets are most likely to be observed in the lower left 
of this figure, and hence may experience significant tidal heating. The recently-announced 
~ 10 Me planet CoRoT-7 b orbits a ~ 0.75 Ms star and hence has a large geometric transit 
probability, but the transit depth was very low, requiring space-based observations to detect 
it. 

Lower mass planets undergo less tidal heating and block less light, i.e. the yellow strips 
in Fig. [I] move to the left, and the horizontal red lines move down. The geometric transit 
probability scales as the sum of the stellar and planetary radii (which scales with mass). 
But since terrestrial planet radii are much smaller than stellar radii, the diagonal red lines 
change negligibly with planetary mass. The IHZ is independent of the planetary mass and 
radius (as long as the planet is massive enough to retain an atmosphere). 

Next we consider the ITHZ for a range of M*, M p , a and e. Figure [2] shows the two HZ 
boundaries for 1 and 10 Me planets (columns) orbiting 0.15 and 0.25 Ms stars (rows) as a 
function of a and e. The color scheme is the same as in Fig. [TJ 

As stellar mass increases, both the IHZ and THZ move to larger a, and the two HZs move 
more or less together. However as planetary mass increases, the THZ may sweep through 
the IHZ, everything else held constant. The IHZ moves to larger a at larger e because the 
orbit-averaged flux increases with e (Williams & Pollard 2002). The tidal heating limits 
move to larger a with larger e because larger eccentricity leads to more heating. 

We can apply these concepts to GJ 581 d (Mayor et al. 2009), whose mass and orbit were 
recently revised to M p = 7.09 Me, a = 0.2184 AU, and e = 0.38. These parameters predict 
the planet is in the Selsis et al. (2007) IHZ. However, using a stellar mass of 0.31 Ms (Udry 
et al. 2007), the tidal heat flux is 0.01 W m -2 = 0.25h min . Therefore our analysis excludes 
this planet from the ITHZ. However, if the planet has at least 0.03 W m -2 of radiogenic 
heating, it could still be habitable. 
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4. Discussion and Conclusions 

We have shown that consideration of tidal heating and stellar flux are crucial when 
assessing planetary habitability. Tidal heating and detectability scale with stellar and plan- 
etary masses, but decrease with increasing a. Hence many planets found in the IHZ may in 
fact be uninhabitable due to intense global volcanism and rapid resurfacing. Consideration 
of tidal heating provides an important refinement in the definition of an HZ. 

We have applied our model to the planet GJ 581 d (Mayor et al. 2009). This planet, 
although in the IHZ, does not reside in the THZ. However, this planet may still possess plate 
tectonics through radiogenic heating. At this point we cannot conclude that GJ 581 d is 
habitable. We encourage future work that examines the radiogenic heating of this planet. 

We have assumed that tidal heating can drive plate tectonics, including subduction. 
For this possibility to be true, heating has to occur deep enough inside the planet to drive 
mantle convection. No models or observations of tidally-driven subduction have been made. 
Such geophysical modeling would be illuminating and could necessitate a revision of the 
ideas presented here. 

Our choice of h m i n , though rooted in observations in the Solar System, could be inap- 
propriate in some cases. The process of plate tectonics is studied by a large and active group 
of geophysicists, and as new results come to light, we may need to change our working defini- 
tion of h m i n . Similarly, our choice for h max may require revision. Planets do not recycle their 
lithosphere uniformly; "cool spots" could persist for long times and be oases for biogenesis. 
Furthermore, planets with unusual compositions or surface features may be able to maintain 
a biosphere even with such large heating. Future geophysical and biological modeling could 
elucidate the true value of h max . 

Although we have generally focused on tidal heating, any heating with the potential 
to drive plate tectonics may stabilize climates. Therefore the ITHZ is really a subset of a 
broader "insolation-geophysical HZ" . It may be that very few planets have h < h min because 
radiogenic heating is generally larger than h m i n (but only for a finite amount of time). 
However, the distribution of radiogenic isotopes in the galaxy remains unknown. Although 
the origin of short-lived radionuclides such as 26 Al in the Solar System is not fully understood, 
models suggest that there should exist a dispersion of several orders of magnitude in their 
in protoplanetary disks, and most planetary systems should actually contain less 26 Al than 
the Solar System (Gaidos et al. 2009; see also Gounelle et al. 2009). If abundances of K and 
U (the main sources of internal heat on the Earth) also have a large dispersion, then many 
exoplanets may not support radiogenically-driven plate tectonics. For such planets, tidal 
heating may be the only available energy source (Williams et al. 1997). 
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The orbital history of a planet, as well as where it is now, may be of critical importance 
when assessing habitability. Most planets inside 0.2 AU have orbits nearly circularized by 
tides from earlier eccentric orbits (Jackson et al. 2008c). As e dropped (over billions of 
years), so did a. In general, e decays more quickly than a (planets tend to move vertically 
downward in Fig. |2} so a planet can form interior to the ITHZ on a high-e orbit which 
is then circularized by tides into the ITHZ. If we detect a planet with such a history, it is 
probably not habitable. The reason is that large terrestrial planets, e.g. ~ 10 Me, will tend 
to hold onto their atmospheres, and hence may be dominated by outgassed species. Thus, a 
planet in the ITHZ today could be uninhabitable due to a volatile past. Spectral analyses of 
transiting planets could reveal evidence of extreme volcanism, even though it is not active 
today. We encourage atmospheric modelers to consider the evolution of such a planet. 

As tides modify the orbit, the eccentricity may eventually reach zero, and Eq. (TjQ) pre- 
dicts no heat flux. Without radiogenic heating, the planet would become uninhabitable. 
However, other planets in the system may perturb the orbit and maintain a nonzero eccen- 
tricity (Mardling 2007), such as in the Galilean satellite system. Models of the formation 
of large terrestrial planets in small orbits suggest such planets may be likely (Raymond et 
al. 2008). Therefore, despite the persistent damping of e by tides, some planets may reside 
in the ITHZ for a very long time. 

As discussed in § 1, tidal effects are but one of several exotic issues that could affect the 
habitability of planets orbiting M stars. For example, the stripping of atmospheres due to 
massive stellar winds (ablation) is also likely a function of M*, M p , and a (Zendejas et al. , 
in prep.). In principle, detailed numerical modeling of these phenomena could identify an 
"ablation HZ" , which could be compared to the ITHZ presented here. The identification of 
an "Insolation-Tidal- Ablation HZ" may guide the wide array of resources that are coming 
on line to detect terrestrial mass planets. Other refinements to the limits of habitability 
incorporating other stellar and planetary phenomena, such as obliquity (Williams & Pollard 
2003; Spiegel et al. 2009), the presence of a large Moon (Laskar et al 1993; but see also 
Barnes & O'Brien 2002), or the presence of water on the planet (Kuchner 2004; Raymond 
et al. 2007b), could be just as useful as the consideration of ablation. 
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Fig. 1. — Boundaries of habitability for a 10 Me terrestrial planet, and associated transit 
detectability levels. The blue region represents the IHZ boundaries for a planet with 50% 
cloud cover (the boundaries change negligibly for e < 0.5). The two yellow strips correspond 
to planets with favorable amounts of tidal heating if e = 0.01 (left strip), or e = 0.5 (right 
strip). Horizontal red lines correspond to different transit depths (with the value of each line 
to the right). Diagonal red lines represent geometric transit probability (with the value of 
each line at the top). 
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Fig. 2. — Habitability boundaries for a 1 Me planet (left column) and 10 Me planet (right 
column) orbiting a 0.15 Ms star (top row) and 0.25 Ms star (bottom row). This figure uses 
the same color scheme as Fig. [IJ 



